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We present further results on a class of sums which involve complex powers of the 

^— V . distance to points in a two-dimensional square lattice and trigonometric functions 

pg ' of their angle, supplementing those in a previous paper (McPhedran et al, 2008). 

We give a general expression which permits numerical evaluation of members of the 

class of sums to arbitrary order. We use this to illustrate numerically the properties 

of trajectories along which the real and imaginary parts of the sums are zero, 

f^ ' and we show results for the first two of a particular set of angular sums (denoted 

C(l, 4to; s)) which indicate their density of zeros on the critical line of the complex 

exponent is the same as that for the product (denoted C(0, 1; s)) of the Riemann 

tG . zeta function and the Catalan beta function. We then introduce a function which 

'T^' is the quotient of the angular lattice sums C(l,4?7i;s) with C(0,l;s), and use its 

(^ . properties to prove that C(l, Am; s) obeys the Riemann hypothesis for any m if and 

only if C(0, 1; s) obeys the Riemann hypothesis. We furthermore prove that if the 

Riemann hypothesis holds, then C(l, Am; s) and C(0, 1; s) have the same distribution 

of zeros on the critical line (in a sense made precise in the proof). 

Keywords: Lattice sums, Dirichlet L functions, Riemann hypothesis 

1. Introduction 

This paper adds to results in McPhedran et al (2008) (hereafter referred to as I) 
on the properties of a class of sums over two-dimensional lattices involving trigono- 

t^^ ' metric functions of the angle to points in the lattice, and a complex power 2s of 

their distance from the lattice origin. There, it was shown that certain of these 
angular sums had zeros on the critical line Re(s) — 1/2, but could not have zeros 
in a neighbourhood of it. 

We derive a general expression which is exponentially convergent and permits 

^ . the rapid and accurate evaluation of the angular sums irrespective of the value 

^ ' of the complex parameter s. We demonstrate the high-order convergence of this 

formula by using it to illustrate a limiting formula for a particular set of angular 
lattice sums. We go on to consider the properties of trajectories along which the 
real and imaginary parts of a class of angular sums are zero, and in particular we 
establish accurate approximations for these trajectories when Re(s) lies well outside 
the critical strip < Re(s) < 1. We give preliminary results on the distribution of 
zeros on the critical line Re(s) = 1/2 of two angular sums, which suggest that to 
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leading order they have the same density of zeros as the product of the Riemann 
zeta function and the Catalan beta function. In the last section of the main body 
of the paper, we introduce a new quotient function, and from its properties we 
prove that the Riemann hypothesis for the product of the Riemann zeta function 
and the Catalan beta function is equivalent to that for all members of a particular 
infinite set of angular lattice sums. We also show that when this Riemann hypothesis 
holds, all members of the class of angular lattice sums have the same distribution 
function for zeros on the critical line, and that distribution function also applies 
to the product of the Riemann zeta function and the Catalan beta function. In 
Appendices we comment on the functional equation satisfied by a class of angular 
sums, and the properties which link angular sums of order up to ten. 

The analytic results presented here are supported by numerical results obtained 
using Mathematica 7.0.1. Formal proofs of certain key properties of the angular 
sums and the location of their zeros will be given here, while others will be presented 
in a companion paper. 

There are two principal motivations for the study presented here. The first 
is that the very general expression for the class of angular lattice sums derived in 
Section 2, and their connections with other angular lattice sums shown in Appendix 
A, enables them to be used in physical applications requiring regularization of 
sums over a two dimensional square lattice. The class of summands which can 
be addressed is wide, as it consists of any function which has a Taylor series in 
integer powers of trigonometric functions of the lattice point angle in the plane, 
and any complex power of its distance from the origin. For example, any summand 
of the type often encountered in solid state physics combining a Bloch-type phase 
factor and a function of distance having a Taylor series could be so represented. 
The second is that we show the angular lattice sums to be connected with the 
product of the Riemann zeta function and the Catalan beta function in a very 
natural way- for example, it seems that the densities of their zeros on the critical 
line are the same to leading orders. In this way, these angular sums may provide 
a new way to forge a link between the Riemann hypothesis and its generalization 
to other Dirichlet L functions, as well as providing a wide class of functions, in 
which identified members obey a priori the hypothesis, and others do not. There 
are interesting parallels between this work and that of S. Gonek (2007), although 
we deal with double sums and Gonek with single sums. 

It should be noted that the proofs presented are not of great technical difficulty, 
and yet some of the results are striking. It thus seems that the study of angular 
lattice sums may offer a favourable context in which properties related to the Rie- 
mann hypothesis may be studied. It is the authors' hope that the results given will 
thus stimulate more mathematicians to recommence the investigation of doubly 
periodic (and even multiply periodic) lattice sums- an area which seems to have 
been unduly neglected for a considerable period. 
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2. An absolutely-convergent expression for angular lattice 

sums 

We recall the definition from (I) of two sets of angular lattice sums for the square 
array: 

where 0pi.p2 = arg(pi + ip2), and the prime denotes the exclusion of the point at 
the origin. The sum independent of the angle 0pi,p2 ^as evaluated by Lorenz (1871) 
and Hardy(1920) in terms of the product of Dirichlet L functions: 

C(0, m; s) = 5(0, m; s) = C(0, 1; s) = 4Li(s)L_4(s) = 4C(s)i_4(s). (2.2) 

Here Li{s) is more commonly referred to as the Riemann zeta function, and i_4(s) 
as the Catalan beta function. A useful account of the properties of Dirichlet L 
functions has been given by Zucker & Robertson (1976). 

It is convenient to use a subset of the angular sums (|2.ip as a basis for numerical 
evaluations. We note that the sums C(n, 1; s) are zero if n is odd. We next derive 
the following relationship for the non-zero sums C(2n, 1; s): 

^ , pf ^ ,_ 2^r{s + n-l/2)C{2s-l) 



r(s + n) ^ ^ \pi 



\ P»"^s+„-i/2(27rpiP2), (2.3) 



where K^{z) denotes the modified Bessel function of the second kind, or Macdonald 
function, with order v and argument z. The general form (j2.3p may be derived 
following Kober (1936) in the usual way: a Mellin transform is used to give 

y ^ 2 2^^ = y 'w^^ / t^+"-'e-*(P?+P^)di. (2.4) 

(Pl,P2) (Pl,P2) 

The Poisson summation formula is then used to transform the sum over p2, giving 



V^ / Pi V^ / Pi 



,^ , b?+P2)'+" ,^ , r(s + 7i) 

\Pl;Vl) (P1.P2) 



r'"""''""'^'""""*- '"' 



We then separate the axial contribution, which for n 7^ comes from P2 = alone, 
and use Hobson's integral 

r i^-ie-P'~«/*di = 2 r^") Ks{2^) (2.6) 

on the remaining double sum. This leads directly to (j2.3p . 

It should be noted that the double sum in (j2.3|) is exponentially convergent. 
Indeed, from relation 9.7.2 in Abramowitz and Stegun (1972), the large argument 
approximation for the Macdonald function of order v is 

K.{^) - /^e-^ (2.7) 
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Figure 1. The modulus of C(2m, 1; s) (red) and 2({2s) (green) as a function of 
s = 1/2 + it, for t e [0, 20], with m = 10 (left) and m = 100 (right). 



This means that the double sum starts to converge rapidly as soon as the argument 
2iTpiP2 everywhere exceeds the modulus of the order s + n— 1/2. In practice, accu- 
rate answers are achieved when sums are carried out over pi and p2 from 1 to P, 
where P ^ \s + n — 1/2\/tt (the precise value of P required being fixed by studies 
of the effect of increasing P on the stability of the result). The representation (I2.3p 
and finite combinations of it thus furnish absolutely convergent representations of 
trigonometric sums from the family (|2.ip and close relatives, for any values of s 
with finite modulus. These representations are easily represented numerically in 
any computational system incorporating routines for the Riemann zeta function of 
complex argument, and Macdonald functions of complex order and real argument. 
(In practice, the Macdonald function evaluations are most time-expensive in the re- 
gion of (pi,p2) values where piP2 has comparable magnitude to |s + n — l/2|/(27r). 
Thus, it is efficient to create a table of these values for pi varying with p2 = 1, run- 
ning up to an argument where the Macdonald function is less than an appropriate 
tolerance times its value for pi = 1.) 

As an example of the numerical efficacy of (12. ip . we consider its use in illustrating 
a limiting property of the sums C{2m, 1; s).We have 



C(2m,l;s) 



E' 



{pI+pIY 



(2.8) 



and as TO — >■ oo we require | cos Op^ ^p^\ — \ for a contribution, i.e. 

lim C(2m, l;s) = 2C(2s). 



(2.9) 



The relationship (j2.9p is illustrated numerically in Fig. [T] For the right-hand side of 
(|2.9p to be accurate, the required order m increases with t, although convergence 
is also slow for t near zero. 

Another angular sum of great importance in this paper can easily be expanded 
in terms of the C(2n, 1; s): 



C(l,4m;s)= ^ 



Pl,P2 



/ COs(4to6>p,,pJ \ Q(n n n 

/ 2 I 2\s = ^^(2, 2to; s) - S{2, 2m; s), 



(2.10) 
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or, in terms of the Chebyshev polynomial of the first kind (Abramowitz & Stegun 
(1972), Chapter 22), 



C(l,4m;s)= ^ 



Pl,P2 



I J 4m (cos C7pi,p2 j 

ipl+plY 



(2.11) 



As the coefiicients of this Chebyshev polynomial are explicitly known, the represen- 
tation (j2.1ip enables any sum C(l, 4to; s) to be expressed as a linear combination 
of sums C{2n, 1; s) with < n < 2m. 

The angular sums just introduced may also be expressed as sums along rays in 
the square lattice. This means that the particular independent values of Op-^ ^p^ oc- 
curring in the lattice are identified, and sums along those directions are performed, 
giving a multiplying factor of C(2s)- The sums along rays may be reduced to those 
occurring in the first octant, which we define as follows: 



Ci = {pi,P2\pi = l,...,oo,l <P2 <Pi,P2 -Lpi}, 



(2.12) 



where the notation p2 -L Pi means the integers are relatively prime. The axis P2 ^ 
and the line pi = p2 are dealt with separately, giving 



C(0,l;s) = 4C(2s) 



l + :i + 2E 



cos 9„ 



pf 



(2.13) 



The sum in (|2.13l) converges down to 3fi(s) = 1, where it diverges, as is to be 
expected from (|2.2p . Other expressions of this type are: 



C(l,4m;s) = 4C(2s) 

C(2,2to;s) = 4C(2s) 1 
and 

5(2, 2m; s) = 4C(2s) 



1 



cos(TO7r) 



2E 



cos(4TO0pj^P2) cos^'* ^Pi,P2 



Oi 



p\' 



cos^(m7r/2) 



2E 

Oi 



cos^ {2m9p^ ^P2 ) cos^'* 9p^ ^p^ 



pf 



sir? {rmi/ 2) 



2E 



sin^(2TO6'p^,P2) cos^" ^pi,P2 



pf 



(2.14) 
(2.15) 

(2.16) 



The last of these formula we give may be compared with 



C(2m, 1; s) = 5(2to, 1; s) = 2C(2s) 



1 



OS + '/Tl 



- + 2E' 



(cos^'"(0p,^P2) + sin^"(0p,,p2)) cos 



^Pl,P2 



Oi 



pf 



There is a sum rule connecting the C(l, 4m; s): 

oo ^ 

^C(l,4m;s) = 7rC(2s)--C(0,l;s). 



(2.17) 



(2.18) 



771 — 1 



This follows from the basic form of the Poisson summation formula (relating a sum 
of cosine functions to a sum of delta functions) . 

The connections between various angular lattice sums grouped in systems with 
order up to 10 are explored in Appendix A. 
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3. Some properties of trigonometric lattice sums 

The functional equation is known (see McPhedran et al. (2004), eqs. 32 and 59) for 
C(l,4m;s): 

G4™(s) = C(l, 4m; ,)£(l±M ^ ^4^(1 - s). (3.1) 

This equation also holds for tti = 0, where it gives the functional equation for 
the product C,{s)L-i{s). It is in fact the m dependence of the functional equation 
(|3.ip which enables the derivation of many of the results in (I) and the present 
paper. As the derivation in McPhedran et al (2004) uses different notation to that 
in subsequent papers and here, we give a brief discussion of the argument leading 
to p.ip in Appendix B. 

This m dependence in (j3.ll) is represented in two related functions: 

r(i-s + 2m)r(s) . , ^^ , , 

^2m{s) = ^^ .^, '; = exp 2102™ s , 3.2 

1(1 — s)l (s + Im) 

where 4>2m{s) is in general complex. Note that J-2m{s) is the ratio of two polynomials 
of degree 2m, with one obtained from the other by replacing s by 1 — s: 

{2m-s){2m-l-s)...{l-s) 
^2m{s) - ^^^ _^^_ ^^jp^ _ 1 _ (1 ^ 5)] . . . [1 - (1 - s)] ■ ^ ' 

We then introduce two sets of rescaled lattice sums: 

C(2, 2m; s) = — J^kL==[C(0, 1; s) + C(l, 4m; s)], (3.4) 

and 

S{2,2m-s) = —-^j^===[C{Q,l-s)~C{lAm-s)\. (3.5) 

Note that this definition means C and S have branch cuts where J-2m{s) is real and 
negative. For example, for J-2{s), the branch cut includes a circle in the (cr, i) plane 
(s = cr + ii), with centre (1/2, 0) and radius a/3/2. 

An analytic function which combines lattice sums is defined by 

A3(2, 2m; s) = C(2, 2m; sf - 5(2, 2m; sf = ^5-^^!— ^C(0, 1; s)C(l, 4m; s). (3.6) 

7r^''-A2m(s) 

It obeys the functional equation 

A3(2, 2m; s) = J-2„(l - s)A3(2, 2m; 1 - s), (3.7) 

and on the critical line s = 1/2 + it, 

A3(2, 2m; s) = [1 - itan((/.2m(s))][|C(2, 2m; s)p - |5(2, 2m; s)W (3.8) 

and 

Im[C(2, 2m; 1/2 + it)] = Im[5(2, 2m; 1/2 + it)]. (3.9) 
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We note from (P?5)) that 

Im A3(2, 2m; - + it) = - taii{<j)2m,ct) Re A3(2, 2m; - + it), (3.10) 

using the notation 02m(l/2 + it) = 4'2m.c{t), a real-valued function. We can take 
the derivative with respect to t of p.lOp . to obtain: 

— ImA3(2,2m;-+it)+tan<?;.2,„,c(i)^ReA3(2,2m;-+it) 

= -;;;J^^'^^^^^3(2,2m;i+if). (3.11) 

COS^(02m,c(i)) 2 

We use the Cauchy-Riemann equations in p. lip , to obtain 

d 1 d 1 

— ReA3(2,2m;-+i<)+tan02m,c(i)iTrReA3(2,2m;-+it) 
oa 2 at 2 

^2m,M .ReA3(2,2m;i+it), (3.12) 



and 



COs2((/)2„,c(0) ' ' 2 



d I d 1 

— Im A3(2, 2m; - + it) - tan ,/,2,„,,(t)— Im A3(2, 2m; - + it) 

Ot Z 0<J I 



*2m,c(i) .ReA3(2,2m;i+it). (3.13) 



COs2(02m,c(i)) 

The equation (I3.12p indicates that, at points where contours of Re A3(2, 2m,;cr + 
it) = intersect the critical line, their tangent vector is given by (1, tan02m,c(O)i 
provided 9A3(2,2m; i + it)/dt ^ 0. (This requirement is that the left-hand side 
of (13.121) can be interpreted as the scalar product of the tangent vector and the 
gradient of Re A3, with the latter having a well-defined direction.) The equation 
p.l3p indicates that the tangent vectors for the contours of Im A3 = at the critical 
line are given by (— tan(/)2m,c(i), l),i.e. they are at right angles to those for the real 
part. 

We note that, for \s\ >> 1, or, more strictly, for \s\ >> 4m^, 

4m^ Sm"* 2m'^i im'^{8m'^ — 1) 

^'"'^'^^^'7^Tj2+is~l/2)^^ ^^"^')"7^T7^+ 6(.- 1/2)3 ■ (3-14) 

Thus, for \t\ >> 1, (I3.12P takes the approximate form 

d 1 2m^ d 1 2m? 1 

— Re A3 (2, 2m; - + it) + ^ _ Re A3(2, 2m; - + it) = -^ Re A3(2, 2m; - + it). 

(3.15) 
This shows that, as t increases, the contours of Re A3(2, 2m; ct -I- it) = strike the 
critical line at ever flatter angles (although the angle increases as m increases). The 
direction of the intersection with the critical line is unique for every point where 
i9A3(2,2m; j+it)/dt 7^ 0. Similar remarks apply to Im A3 (2, 2m; a -I- it) = 0, where 
the equivalue contours cut the critical line at a direction given by the tangent vector 
(— 2?7t^/t, 1), and thus their gradient at the point of intersection increases with t. 
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Figure 2. Null contours of the real part (red) and imaginary part (blue) of A3 (2, 2; cr + ii), 
with (left)(T £ [—0.5, 1.5], and t G [0, 10]. On the right we show part of the left panel, with 
now t £ [0.0,2.0], to give an undistorted geometric image. On the green circle, J-2{s) is 
real and negative. 

4. Equivalue contours of Re A3 and Im A3 

In Figs. [2]|l]we show contours on which the real and imaginary parts of A3(2, 2; s) 
are zero. Fig. [5] gives some detail of the region near a = 0.5 for i G [0.1, 10], and a 
further more detailed region on an undistorted geometric scale (including the semi- 
circle on which J^2is) is real). Fig. |3] gives a more global view of null contours for t 
ranging up to 20. Fig. |4] gives the detail of the null contours for t near two values at 
which contours nearly touch. We introduce the notation 14 for a zero on a = 1/2 of 
C(l, 4; s), +1 for a zero of ({s) and —4 for a zero of L_4(s). Then the zeros evident 
in Figl3]are categorized as: -4, 14, 14, -4, 14, -4, +1, 14, -4, 14, -4, 14, 14. 

The null contours have a number of interesting properties. Firstly, we can see 
that contours for the real part do indeed strike the critical line at a small angle 
to the a axis, which decreases as t increases, while the contours for the imaginary 
part strike the critical line almost vertically. Secondly, the contours for the real and 
imaginary parts intersect the critical line simultaneously, except for one point. This 
is at (1/2, •\/3/2), where, as remarked in (I), 02(1/2 + 1^3/2) = 7r/2, which permits 
ReA3(2,2;s) to be zero, while ImA3(2,2;s) is non-zero (see (|3.8p ). 

We can see four null contours for the real part and five null contours for the 
imaginary part intersecting the axis t = 0. If we use the expression 



A3(2,2;s) 



T{s)s{s + l) 
(2-s)(l-s)7r2 



:[C^{2,2-s)-S\2,2;s)], 



(4.1) 



we find that A3 (2, 2; s) has a second order pole at s = 1, and first order poles at 
s — and s — 2. Numerical investigations show that near these points 



A3i2,2;l + S) 



-1.59643 



(4.2) 
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Figure 3. Null contours of the real part (red) and imaginary part (blue) of 
A3(2,2;cr + it), with a G [-4.5,11.5], and t G [0.1,20]. 



and 



^3(2,2:6) 



-0.798212 



A3(2,2;2 + <5) 



1.16981 



(4.3) 



We note that ReA3(2,2;s) has zeros at the two first-order poles, where the tra- 
jectories approach the poles broadside (i.e., parallel to the t axis), and two null 
trajectories approaching the second-order pole at 45° and 135° to the a axis. 
ImA3(2,2;s) has a zero at the second-order pole, where the null trajectory ap- 
proaches broadside. The other null trajectory crossings of the real axis occur at 
0.29782, 1.67735, -2.65568 and 4.21422. The last two of these are associated with 
minima of Re A3 (2, 2; a). 

We show in Fig. [5] the behaviour of the null contours for a large enough (ct > 4) 
to enable us to accurately approximate trigonometric sums by their first few terms: 



C(2,2;s)2-5(2,2;s) 



16(1 



36 , 

r2.s+2 ' 



(4.4) 
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Figure 4. Detail of the null contours of the real part (red) and imaginary part (blue) of 
A3 (2, 2; o- + if), for t near 10 (left) and near 18.5 (right). 




Figure 5. Equivalue contours of the real part (red) and imaginary part (blue) of the form 
gU for A3(2, 2; cr + ii), with a G [3.5, 7.5], and t e [0.1, 10]. 



SO that 



A3(2,2;s) 



i6r(i-s)r(s + 2)r(s) 



r(3-s)7r2 



1 36 , , , 

(l + TI + T^JTl)- (4-5) 



For a beyond 4.21422, we see from Fig. [5] that neither the null trajectories of the 
real part or the imaginary part of A3(2,2;s) can attain the real axis. In Fig. |6l 
we show for comparison the null trajectories associated with the prefactor term in 
square brackets in (|4.5p . It is evident from a comparison of Figs. [S] and [S] that the 
null trajectories are given accurately by the prefactor in the range of cr shown. Using 
Stirling's formula (Abramowitz and Stegun (1972)), we can place the prefactor in 
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Figure 6. Null contours of the real part (red) and imaginary part (blue) of the prefactor 
term in the expression (P3)) for A3(2, 2;cr + ii), with a £ [3.5, 7.5], and t e [0.1, 10]. 

exponential form for jsj large, and find the constraint for null trajectories to be 



Im[2slog(-) -2s 



log(s) + 



25 
6s 



{n -\ — )tt, or WK, 



(4.6) 



where on the right-hand side the first value is for real-part zeros, and the second 
for imaginary-part zeros {n being an integer). From (I4.6|) we find that the equiv- 
alue contours alternate for a large, with null contours for the real part sandwiched 
between those for the imaginary part, and vice versa (with each trajectory corre- 
sponding to a first-order zero). Each contour tends to zero as 1/ log(CT/7r) as a — > oo. 



Note that we may use the equation (j3.7p to deduce relations governing the phase 
of A3(2, 2to; s). Dividing (13.71) by its conjugated form, we obtain 



arg(A3(2, 2m; s)) + arg(A3(2, 2m; I - s)) = - &ig[F2^{s)). 



(4.7) 



This agrees with (13. 8p and p.lO[) when Re(s) = 1/2. Furthermore, on the right-hand 
side of Fig. [3]we know that arg(A3(2, 2; s) increases monotonically with i, from its 
value of zero on the real axis. This enables us to assign values of the constant phase 
of the null contours there: 7r/2, tt, — 7r/2, 0, 7r/2, etc. Using (|4.7p and (|3.14p . we see 
that the monotonic increase of arg(A3(2, 2; s) with t on the right in Fig. [S]) forces 
a monotonic decrease of arg(A3(2, 2; s) with t on the left. We can then assign the 
phase values of the null contours there, again starting from zero on the real axis: 
-7r/2, -TF, 7r/2, 0, -7r/2, etc. 

We can use this information to understand the behaviour of the null contours 
shown in Fig. [3l As the null contours of Re(A3(2, 2m; s)) pass through a simple 
zero, their phase must change by tf, from 7r/2 to — 7r/2, or vice versa. However, we 
note that the relevant contours on the left in Fig. |3] which are almost symmetric 
to those on the right have opposite phase values. This means the null contours of 
Re(A3(2, 2m; s)) vary smoothly as they cross the critical line. For the null contours 
of Im(A3(2, 2m; s)) the situation is different: the phase change of tt forces them 
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to avoid their almost symmetric counterpart on the left and "jump up" a contour 
as they pass through a zero on the critical line. These remarks are in accord with 
the derivative estimates at the critical line (see equation (j3.15p and subsequent 
discussion). 

5. Some properties of zeros 

The properties of the angular sums given in the preceding sections enable simple 
proofs to be given of some properties concerning the location of zeros of angular 
lattice sums. These properties form an interesting counterpoint to those of the 
zeros of Epstein zeta functions (see Bogomolny and Lebouef (1994) and references 
[26-30] of that paper) . The Epstein zeta functions are characterized as having an 
infinite number of zeros lying on the critical line, but with many zeros lying off 
that line and with almost all zeros lying on the critical line or in its immediate 
neighbourhood. We have already seen in the example of Section 2 that this third 
property of Epstein zeta function zeros is not shared by those of the trigonometric 
double sum C(2m, 1; s) for large m. 

Theorem 5.1. The trigonometric sums C{2,2m; s) and S{2,2m; s) have no zeros 
for s on the critical line in the asymptotic region t >> Am? which are not zeros of 
both, and these will then he simultaneous zeros of C{0, 1; s) and C(1,4to; s). 



Proof. We combine equations (13. 4p and (13. 5|) . to obtain 

(5.1) 



C(2,2m;s) + 5(2,2m;s) "^ 



^CiO,l;s) 



As the term in square brackets is real on the critical line, from the functional 
equation (j3.ip . we see that 



arg[C(2, 2m; 1^+^) + S{2, 2m; - + it)] = -^.n-dt) + 



(5.2) 



The second term on the right-hand side of (15. 2p is chosen according to the sign of 
the term in square brackets in (j5.ip . We can also use (j3.9p to deduce that 



arg[C(2, 2m; \ + it) - S{2, 2m; - + tt)] = 



(5.3) 



Note that the term on the right-hand side in (15.31) incorporates the phase of the 
real function on the left-hand side, and is not necessarily the same as the second 
term on the right-hand side of (J5.2I) . 

If we suppose that 5(2, 2m; ^ + it) = 0, wc require from (j5.2p and (|5.3p that 
4>2m:c{t) = n-K for somc integer n. We know from (I3.14p that this cannot occur. 
Exactly the same argument applies if C(2, 2m; ^+ it) — 0, so neither function can 
be separately zero on the critical line in the asymptotic region. If both are zero 
then, trivially, C(0, 1; s) = and C(l, 4m; s) = 0. D 

We now consider the phases of S{2, 2m; ^ + it) and C(2, 2m; ^+it), particularly 
in the neighbourhood of zeros C(0, l;s) or C(l,4m,;s). For brevity, we will adopt 
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the following notations: 

\Ci2,2m;^+2t)\ = \Cl argC(2, 2m; 1 + z<) = 9^, (5.4) 

and 

\Si2,2m;^+it)\^\S\, arg^(2, 2m; i + z<) = 9,. (5.5) 

Then on the critical line, we have for C and S that 

C{2,2m;-+it) = |C|(cos9c + isin9c), S{2,2m; - + it) = |S'|(cos9s + isin9s), 

(5.6) 
and ^ ^ 

|C|sin9c= |S'|sin9s. (5.7) 

The relationships between these phase angles are exemplified in Fig. [3 where the 
variations with t of the functions cot(9c) and cot(9s) are shown for the case m = 1. 

Theorem 5.2. The phase angles 9c and 9s obey the following equations on the 
critical line 

COt(9c) + C0t(9s) = -2 COt(02m,c(t), (5.8) 

and 

I cot(9.) - cot(9.)| = 1^,,^'^S!;^";''^' ,., ■ (5.9) 

\C[0,l; s)\\sm(4>2m,S)\ 

Proof. We prove these relationships using the diagram of Fig. [51 in which complex 
quantities are represented by vectors. For the case shown, the angle 9^ runs from 
the extended line OA round to AD, while 9c runs from line OA round to OM, with 
the former angle exceeding the latter. The vectors representing C — S and C + S 
are respectively OA and OD, while the angle between OA and OD is tt — (j)2m,c{t)- 
We apply the sine rule to the triangle OAB, to give 



2|C| _ \C' + S\ _ \C~S\ 



sin(02m,c(i)) sin(9c) sm{Qc + 4>2m.ci't))' 
We next apply the sine rule to the triangle OAD, to give 

2|^| _\C + S\_ \C-S\ 



sin(02m,c(O) sin(9s) sin(9c + 02m,c(O - tt) ' 



(5.10) 



(5.11) 



Using the first two elements of (j5.10p and (|5.1ip we recover the result that |C|/|S'| = 
sin(9s)/sin(9c). Using the second and third elements, we find 

\C — S\ 

cos(02m,c(O) + cot(9c) sin((/)2™,c(t)) = -p- -j, (5.12) 

\C - S\ 



and 



\C — S\ 

cos(02m,c(i)) - C0t(9s) sin((/)2m,c(i)) = -^ -" ■ (5.13) 

\c — s\ 
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Subtracting (I5.12p and (|5.13p . we obtain (|5.8I) . while adding them and using 

\C-S\ _ \C{l,4m;s)\ 
\C+S\ \CiO,l;s)\ 



(5.14) 



we obtain (|5.9p for the particular case of Qg > ©c and sin(02m,c(i)) > 0. The form 
given for (J5.9I) is true whatever the inequality between Qg and Qc, and whatever 
the sign of sm{(j)2,n,c{t))- □ 

The curves of Fig. [7] giving cot{Qc) and cot(9s) on the critical line as a function 
of i for m = 1 illustrate the relationship (15.81) . One may be obtained from the other 
by reflection in the line giving the variation of — cot{(j)2m,c{t)- They intersect when 
©c = ©s, at a zero of C(l,4m; ^ + it). The curves of cot(©c) and cot(9s) also 
go off to infinity at zeros of C(0, 1; ^ + it), with their movement being in opposite 
directions in accord with (j5.8p . In Fig. [7] we see examples of there being none, one 
or two zeros C(l, 4m; | + it) between successive zeros of C(0, 1; 5 + it). 

Using (|5.6p and (|5.7p . we can rewrite (|5.2p as 

SsinGcSinGs ,, , , 

^— —— — — = -tan(/.2m,c i • (5.15) 

sm(ec + 6s) 

With a prime denoting t derivatives, we obtain by differentiating (j5.15p 

(2 cos Be sin 6^)6^ + (2 sin 9^ cos 9^)91 = 
02m,cW .sin{Q, + ec)-tan<l>2m.At)cos{Qs + Qc){e's + Qc)- (5.16) 



COS'^ 4i2-mAt) 



We now consider the values of 9^, 9c and their t derivatives in the neighbour- 
hood of zeros of C(1,47tt,; ^ + it) or C(0, 1; ^ + it), which are not zeros of both 
functions. Close to a zero of C(l, 4m; ^ + it), we require 

|C|cos9c-|5|cos9s^0, |C| ~> |5|, 9^ ^ 9^, (5.17) 

where \C\ 7^ 0, \S\ 7^ at the zero by Theorem [ET] Using the last of ((^Tf)) in ((??T5|) 
and (|5.16p . we find at the zero of C(l, 4m; i + it) that 



9s = 9c = -4>2niAt) + 



9^ + 9; = -2<j>'2„,At)- (5-18) 



Close to a zero of C(0, 1; | + it), we require 
|C|cos9c + |S'|cos9s-^0, 1(71 sin9c +15^1 sin 9s ^0, |C| ^ IS']. (5.19) 

From (|5.19p . at the zero 

9; = -9^. (5.20) 



9c = 7r + 9s = 



The next result relates to the independence of the functions A3 (2, 2m; s) on the 
critical line, again in the asymptotic region. 
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Figure 7. Variation of the quantities — cot((^2,c(i)) (red), cot(0c(i)) (blue) and cot(0s(t)) 
(green) with t £ [9.5,20] (left) and [20,30.5] (right), where the angles 0c and 0s pertain 
toC?(2,2;s) and S(2, 2; s) . 
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Figure 8. The complex quantities (7(2, 2; s) (OM) and 5'(2, 2; s) {AM) represented as 
vectors, together with their sum (OD) and difference (OA), with the last two respectively 
proportional to C(0, 1; s) and C(l, 4; s). 

Theorem 5.3. // .s is on the critical line in the asymptotic region t >> AP and 
ai, (3 are non-zero reals, then 



aiC{l, 41; s) + (3C{0, 1; s) ^ 0, unless C(l, 41; s) = = C(0, 1; s) 



(5.21) 



// s is on the critical line in the asymptotic region t >> 4niax(Z^, tti^) and ai, am 
are non-zero reals with I ^ m, then 

a;C(l, 4Z; s) + a,„C(l, 4rn; s) 7^ unless C(l, 4/; s) = == C(l, 4m; s). (5.22) 



Proof. We suppose 

aiC(l, 4/; s) + /3C(0, 1; s) = with C(0, 1; s) 7^ 0. 
Then 



6(0, l;s) 



g(l,4/;5) 
C(0,l;s) 



(5.23) 
(5.24) 
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for m an integer. 

We now use the functional equations (|3.1[) for C(0, 1; s) and C{l,Am; s), which 
on the critical line with 1 — ,s = s enable their arguments to be deduced: 

argC(0, l;s) = -[argvr'* + argr(s)] , argC(l,4r7i; s) = -[argTr" + argr(s + 2to)]. 

(5.25) 

Using (j5.25|) we arrive at the estimate for s on the critical line in the asymptotic 

region 

2m^ 1 

arg C(0, 1; s) - arg C(l, 4m; s) = ttitt + 0( — ). (5.26) 

t t^ 

Using (Hiini), (|?^ requires 

2^2 1 

— + O(^) - /tt = mTT, (5.27) 

which is not possible. This proves (|5.2ip . 
Consider next (|5.22p . and suppose 

a/C(l, 4Z; s) + a^C(l, 4m; s) = with C(l, 4m; s) ^ 0. (5.28) 

Then this requires 



arg 



_C(l,4m;s) 
for an integer p. However, from (I5.26|) . 



PTT, (5.29) 



argC(l, 4?; s) - argC(l, 4m; s) ^ ira - 1)ti + -^ -^ + O(-). (5.30) 

Hence (j5.28p requires 

^—^ l + 0{^)^[p + l-m)ir, (5.31) 

which is not possible if ^ 7^ m. This proves (j5.22p . D 

Corollary 5.4. The sums C(2, 2, s), 5(2, 2, s) and C(4, 1; s) /laue no zeros s on the 
critical line in the asymptotic region t >> 4 which are not zeros of both C(0, l;s) 
and C(l, 4; s). 

Proof. We apply TheoremlOland equation ((OT|) . Now the sums C(2, 2, s), 5(2, 2, s) 
and C(4, 1; s) belong to the system of order 4 discussed in Appendix A, which is gen- 
erated by linear combinations of C(0, 1; s) and C(l, 4; s). All such combinations are 
of the form (j5.2ip . and are thus zero only if both C(0, l;s) = andC(l,4;s) = 0. D 

We can also simply establish a corresponding result for the derivative A3 (2, 2m; s) 
for s on the critical line. 

Theorem 5.5. The derivative function A3 (2, 2m.; s) has no zeros for s — 1/2 + it 
in the asymptotic region t >> 4m^ which are not zeros 0/ A3(2, 2m;s) of multiple 
order. 
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Figure 9. The real part (red), imaginary part (blue) and modulus (green) of 
A3(2, 2; 1/2 + it), with t G [0.1, 20] (left) and [20, 40] (right). 

Proof. We recall the result ((3^ for s = 1/2 + ii: 

A3(2, 2m; 1/2 + it) = [l-itan(02m,c(i))][|C(2, 2m; l/2 + itf~\S{2, 2m; l/2 + it)\^], 

(5.32) 
where (/)2m,c(i)) is a real-valued function. Taking the derivative of A3(2, 2m.; 1/2+it) 
with respect to 1/2 + it and separating into real and imaginary parts, we find 



A3(2,2m;l/2 + ii) = - 



COs2 02m,c(i) 



[|C(2, 2m; 1/2 + it)\^ - \S{2, 2m; 1/2 + it)\^' 



+ tan02m,c(i)i-[|C(2, 2m; 1/2 + it)\'^ - [5(2, 2m; 1/2 + it)\^] 
at 



^i — [|C(2, 2m; 1/2 + it)\'^ - [5(2, 2m; 1/2 + it)\'^]. (5.33) 



For A3(2, 2m.; 1/2 + it) to be zero, we find from the imaginary part of (j5.33p that 
d 



dt 



[\C{2, 2m; 1/2 + it)\^ - \S{2, 2m; 1/2 + it)\' 



0. 



(5.34) 



We use (|5.34p in the real part of (J5.33L and, assuming ta.n(f>2m,c{t) is non-singular 
and 02^ ^{t) 7^ 0, we complement (J5.34I) with 



[|C(2, 2m; 1/2 + it)\^ - \S{2, 2m; 1/2 -f- it)\'^ 



0. 



(5.35) 



Since both the assumptions we have just mentioned are true in the asymptotic 
region, we have proved the theorem. D 

We provide plots of A3 (2, 2; 1/2 -I- it) in Fig. [3] These show that the derivative 
is non-zero in the non-asymptotic region as well as in the asymptotic region. The 
numerical value of AJ5(2, 2; 1/2) is 0.918604. 

We next consider the properties of trajectories of constant phase which follow 
from the assumption that the Riemann Hypothesis holds for A3 (2, 2m; s). As can 
be seen from Fig. [TUl lines of constant phase between zeros of A3(2,2;s) do not 
in general cross the critical line, but asymptote towards it, and their configuration 
is arranged about a zero of the derivative 9|A3(2, 2; cr -|- it)\/dt or equivalently of 
(9argA3(2,2;cr + it)/acr. 
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Figure 10. Contours of constant phase for A3(2, 2;s) in the region around its zero at 
s = 0.5 + il9.80599, with a second zero evident at s = 0.5 + i21. 02204. The plot at right 
shows detail around the zero of the derivative 9 log |A3(2, 2; cr + it)\/dt. 

Theorem 5.6. Given the Riemann hypothesis holds for A3{2,2m; s), then lines 
of constant phase coming from a = — oo can only cut the critical line at a zero of 
A3(2,27Ti;s) or o/ (9arg A3(2, 2; cr + it)/9(T. Those lines of constant phase coming 
from a = oo cutting the critical line at a point which is not a zero of A3 {2, 2m; s) or 
o/9arg A3(2, 2; a + it)/da must curve back and pass through a zero 0/ A3(2, 2m; s). 



Proof. The assumption of the Riemann hypothesis holding enables us to say that 
lines of constant phase coming from a = —00 do not intersect before the critical 
line. Their phase monotonically decreases as t increases in cr < 1/2, while the 
phase of A(2, 2m; 1/2 + it) monotonically increases with t. Thus, groups of constant 
phase lines coming in from cr = — 00 cannot cut the critical line, except at a zero 
of A(2, 2m; s). Isolated trajectories passing through points where 9arg A(2, 2; cr + 
it) /da = are allowed. Such special trajectories separate lines of constant phase 
which curve up as cr — > 1/2 and they move towards a zero of A(2,2m;s), from 
trajectories which curve down as cr —>• 1/2. 

Now consider sets of lines of constant phase approaching cr = 1/2 from the right. 
Their phase increases with i, but if they were able to cross the line cr = 1/2 and 
progress towards cr = —c)o it would have to decrease with t. Thus they must turn and 
run alongside a — 1/2, with lines above the special trajectory curving up towards 
a zero of A(2, 2m; s), and those below curving down towards a zero. Alternatively, 
they can cut the critical line at a point which is not a zero of A3 (2, 2m; s) or of 
9 argA(2, 2; cr + it) /da. Such trajectories supply the lines of constant argument 
required for generic points on the critical line, and must return back to cut the 
critical line at a zero of A(2, 2m; s). The region where they cross into a < 1/2 is 
bounded on the left by a line of constant phase connecting a zero of A3 (2, 2m; s) 
with an adjacent zero of d arg A(2, 2; cr + it) /da. D 

Note that the region where zeros of 9 arg A (2, 2; cr + it) /da exist is confined 
to the neighbourhood of the critical line, since this partial derivative being zero 
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corresponds to a horizontal segment on a line of constant phase (see, for example, 
Fig. [3]). However, we may use the line of constant phase passing through the point 
on the critical line where i9arg A(2, 2;a + it) /da = as the separator between lines 
of constant phase going to the zero above this line from those going to the zero 
below (given the Riemann hypothesis is assumed to hold). 

Theorem 5.7. Given the Riemann hypothesis holds for /S.z{2,2m\s), then there 
exists one and only one zero of d avg/S.y,{2,2]a + it)/ da on the critical line between 
two successive zeros of /S.j,{2^2m] s). 

Proof. We consider the analytic function 



log A3(2, 2m; s) = log |A3(2, 2m; s)\ + i arg A3(2, 2m; s). 



(5.36) 



The real part of this function goes to — oo at any zero of A3(2,27n;s), and 
increases away from these logarithmic singularities. It must have at least one turning 
point between successive zeros. By the Cauchy- Riemann equations, such a turning 
point is a zero of 9arg A3 (2, 2; cr + it) /da. 

Next, suppose there two or more zeros of (3arg A(2, 2;a + it) /da between suc- 
cessive zeros So = 1/2 + ztg and si = 1/2 -f iti. Denote the upper two of such 
derivative zeros by s■^, = 1/2 + it.^, and s** = 1/2 + it^.^,. As we have seen, each of 
these has constant phase lines coming from cr = — 00 and passing through it, around 
which constant phase trajectories reverse their course. Those above s* curve up to 
si as they approach the critical line, while those below it curve down. They cannot 
cross the constant phase line passing through s**, nor can they cross the critical 
line. They must then head back to ct = —00, where they will breach the mono- 
tonic nature of the variation of arg A3 (2, 2m; s) with t. Thus, this situation cannot 
arise. D 

We conclude this section with an investigation of the structure of lines of con- 
stant phase which cut the critical line and turn back to ct = 00 thereafter. We start 
with the expansion of arg A3(2, 2]a + it) around a point s* = cr* -|- it* = 1/2 -|- it,, 
in the asymptotic region where (9arg A3(2, 2; cr + it)/da = 0, which is of the form: 



2m 
argA3(2,2;CT + it) = arg A3(2, 2; 1/2 + ii*) + -^(a - 1/2) 



+- 



2m^ 



t? 



(i-t*)- 



92 
df' 



tl 



log|A3(2,2;(7 + ii)| 



{a — a^,){t — f*) + 



2m^ 
~tf 



{t-t,f 



(5.37) 



Here we have employed the asymptotic estimates for arg A3 (2, 2; 1/2 -I- it) based 
on (j3.14p . The corresponding trajectories of constant phase may be shown to be 
rectangular hyperbolae, with their centre at 



CT-1/2 = 






1 



1 r d 



[^log|A3(2,2;a + zi) 
1 



1' 



1 r_ai 



[^log|A3(2,2;a + .t)|]^^^^ 






m 
2a 



/2m4 



4™" V ** 



df^ 



log|A3(2,2;a + 2f)| 



(5.38) 
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Figure 11. Contours of constant phase for A3(2, 2;s) in the region around its zero at 
s = 0.5 + il9. 80599. The thick line is the contour on which 9arg A3(2, 2; a + it)/dt is zero. 
(a)The fine lines correspond to the phases 3.04,3.04076, 3.0410, 3.0412, 3.0414, 3.0416, 
3.0418, 3.0419, 3.04198, 3.042, 3.0425, 3.043, 3.044, 3.045, 3.046, 3.047 and 3.048. (b) The 
thick dashed hne is the contour on which 9 log |A3(2, 2;a + it)\/dt is zero. The contours 
are -0.1090, -0.1080, -0.1072, -0.1068, -0.1064, -0.1060, -0.1056, -0.1055, -0.1054, -0.1053, 
-0.1052, -0.1048, -0.104, -0.103, -0.102, -0.1015 and -0.10. 



The second derivative factor in (j5.38p is always negative, and tends to be much 
larger than the terms involving powers of 1/i*. Thus, we see that the centre of the 
hyperbolic trajectories of constant phase will always lie to the left of the critical 
line, with its ordinate very close to i*. This displacement of the centre into a < 1/2 
creates the region in which phase lines can cut through a = 1/2 and return to a > 
1/2 via passage through a zero of A3(2, 2; s). Note that the centre corresponds to a 
point at which two lines of constant phase intersect; thus, it must have derivatives 
of phase along two independent lines which are zero. It therefore is a point at 
which both 9arg A3(2, 2; a + it)/dt and 9arg A3(2, 2; a + it) /da are zero. From the 
Cauchy-Riemann equations, it also is a point at which 51og |A3(2, 2; a + it)\/dt and 
91og |A3(2, 2; a + it)\/d(T are zero- i.e., it is a point of extremum for both amplitude 
and phase. 

These hyperbolic centre points are locations at which A3 (2, 2; s)' = 0. The fact 
that they lie to the left of 5ft(s) = 1/2 is is similar to the property proved by Speiser 
(1934) for C(s). 

The curves of constant phase of A3(2, 2; s) given in Fig. [TT] illustrate some more 
of their general characteristics in the neighbourhood of zeros. For Fig. [TT] (a), at 
the top of the figure, we see the centre of the hyperbolic phase curves; through this 
centre passes a curve on which 9argA3(2, 2; cr + it)/dt — 0. This curve connects 
the hyperbolic centre to the zero of A3(2,2;s) below it, and also continues to the 
zero above; at each zero it is tangent to the t axis. The curve marks points where 
the constant phase lines have vertical slope. The lowest line of constant phase on 
the right has vertical slope when it arrives at the zero of A3(2, 2; s); it corresponds 
to 3.04076, which is the value given by (13.81) for the phase on the critical line 
just above the zero. A second important line has the phase 3.04198, which is the 
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phase corresponding to the centre of the hyperbola. This hne again connects the 
hyperbohc centre to the zero. All lines whose phase lies between these values come 
in from the right, cross the zero line of 9arg A3(2, 2; u + it)/dt, and curve back to 
pass through the zero of A3(2, 2; s). Where they lie to the left of the zero line of 
(9arg A3(2, 2; ct + it)/dt, their phase increases as t decreases; where they lie to the 
right, it increases as t increases. Curves coming in from the left all pass through 
the zero without crossing the zero derivative line; their phase always increases as t 
decreases. 

For Fig. [TT] (b), we show the phase contours below the zero, in the region down 
to the next zero. The phase at the centre of the hyperbolae is -0.1055, while the 
phase just above the bottom zero and just below the upper zero are respectively 
-0.107604 and -0.101134. Note the dashed line passing through the centre, along 
which (91og |A3(2, 2;(T -I- it)\/dt = 0. From the Cauchy-Riemann equations, this 
is perpendicular to the solid line defined by i9arg A3 (2, 2; cr -I- it)/dt = (a fact 
disguised by the different scales on the horizontal and vertical axes). 

6. Distributions of zeros 



We return to the left-hand side of expression (j4.6p . in which we replace s hy a + it, 
and expand assuming |t| >> a, with a large enough to ensure accuracy of (|4.6p . 
The result is 

2tlog(i) - 2i(log7r + 1) + 7r(CT - -) + ^(1 - 2a). (6.1) 

As for each increment of tt of this expression we get one null line of the real part 
of A3 and one of the imaginary part, and these intersect at cr = 1/2 to give one 
zero there (assuming the Riemann hypothesis holds for A3(2, 2m; s)), we can divide 
(|6.ip by TT, and regard the result as a distribution function for zeros of A3: 

NAA<^,t) = - log(t) ^-{l + \og7:)+a-l + -(l-2a). (6.2) 

TT TT 2 TTl 

Now, from Titchmarsh and Heath-Brown (1987), the distribution function for 
the zeros of the Riemann zeta function on the critical line is 

Ndl^t) ^^logit)-^il + \og{27:))+0(logt). (6.3) 

We complement this with the numerical estimate from McPhedran et al (2007) for 
the distribution function of the zeros of L_4(s): 

N-a{\, i) = ^ log(i) - ^(1 + log(7r/2)) + 0(log t). (6.4) 

Adding (J6.3D and (|6.4p we obtain the distribution function for the zeros of C(0, 1; s) 
(see dSH)): 

A^co,i(i = - log(t) - -(1 + log(7r)) + O(logi). (6.5) 

Z TT TT 

When we compare this with (j6.2p . and use the equation 

NA:A\.t) ^ NciA{\,t) + Nca,^{]^,t), (6.6) 

Article submitted to Royal Society 



22 



R.C. McPhedran and others 




Figure 12. The distributions of the differences between successive zeros of C(0, 1; s)(left) 
and C(l, 4; s) (right) on s = 1/2 + it, for t < 300. 

it suggests the hypothesis that the distribution function of zeros of C(l,4; s) is the 
same as that of (|6.5I) , to the number of terms quoted: 

NciA(Lt) = Nco.iiLt) = - logW - -(1 + logW) + O(logt). (6.7) 

Strong numerical evidence supporting this is given in Table (TJ which also shows 
zero counts for C(l, 8; s) and C(l, 12; s). Note that the numbers of zeros found for 
C(l,4; s),C(l,8; s) andC(l, 12; s) are virtually the same. This rules out any variation 
with increasing order similar to that of Dirichlet L functions, where increasing order 
results in significant increases in density of zeros (compare (j6.4p and (j6.3p . or the 
second and third columns of Table 1). 

Comparing the data of Table 1 with the discussion in Bogomolny and Lebouef 
(1994), we can see that the split up of N{t) into averaged parts given by expres- 
sions like (J6.31l6.7p and oscillating parts applies to C(0, 1; s) and to the C(l, 4m; s). 
However, it would be value to extend the numerical investigations of Table 1 to 
much higher values of t, to render the characterization of the oscillating term more 
accurate. Such an extension may require the development of an alternative algo- 
rithm to that based on (|2.3p . which will probably become unwieldy for values of t 
of order 10^- 10^ 

In Fig. [1^ we compare the distributions of the differences between zeros on 
the critical line for C(0, 1; s)(left) and C(l, 4; s) (right). The distributions are quite 
different, even with this modest data set. Bogomolny and Lebouef (1994) have 
studied the case of C(0, l;s) using 10,000 zeros after t = 10^, and contrast the 
distribution for C(s)L_4(s) with that for each function separately. The separate 
factors in fact have distributions like that that for C(l, 4; s). The function compared 
with the histogram in the right of Fig. [12] corresponds to the Wigner surmise, 
which (Bogomolny and Lebouef, 1994, Dietz and Zyczkowski, 1991) for the unitary 
ensemble takes the form 



P(S) = ^ exp 



-452 



Here the separation between zeros has been rescaled to have a mean of one. Bo- 
gomolny and Lebouef (1994) comment that the left distribution is that of an un- 
correlated superposition of two unitary ensemble sets. Note as one indicator of this 
that there is not the same pronounced tendency for the probability to go to zero 
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Table 1. Numbers of zeros of ({1/2 + it), L_4(l/2 + ii), C(l,4; 1/2 + it),C(l, 8; 1/2 + it) 
and C(l, 12; 1/2 + it) in successive intervals oft. 

















t 


nc 


n_4 


nci4 


n_4 -1- nci4 


nci8 


ncii2 


0-10 





1 


2 


3 


2 


3 


10-20 


1 


4 


5 


10 


5 


5 


20-30 


2 


5 


6 


13 


7 


7 


30-40 


3 


4 


8 


15 


8 


8 


40-50 


4 


6 


8 


18 


8 


8 


50-60 


3 


5 


9 


17 


9 


9 


60-70 


4 


6 


9 


19 


10 


10 


70-80 


4 


6 


11 


21 


10 


10 


80-90 


4 


7 


11 


22 


11 


10 


90-100 


4 


6 


10 


20 


10 


12 


0-100 


29 


50 


79 


158 


80 


82 


mi,mi,mi 


28 


50 


78 


156 






100-110 


4 


7 


11 


22 


11 


10 


110-120 


5 


7 


12 


24 


12 


12 


120-130 


5 


7 


12 


24 


12 


12 


130-140 


5 


7 


12 


24 


12 


12 


140-150 


4 


7 


12 


23 


12 


11 


150-160 


6 


7 


12 


25 


12 


13 


160-170 


6 


7 


13 


26 


13 


12 


170-180 


6 


8 


13 


27 


13 


13 


180-190 


5 


8 


13 


26 


12 


14 


190-200 


5 


7 


13 


25 


14 


13 


0-200 


80 


122 


202 


404 


203 


204 


(lOf.llOll.lIOl) 


78 


122 


200 


400 






200-210 


6 


8 


13 


27 


12 


13 


210-220 


5 


8 


14 


27 


15 


14 


220-230 


6 


8 


13 


27 


13 


13 


230-240 


6 


8 


14 


28 


14 


14 


240-250 


6 


8 


14 


28 


14 


13 


250-260 


6 


8 


13 


27 


14 


15 


260-270 


6 


8 


15 


29 


14 


13 


270-280 


6 


8 


14 


28 


14 


15 


280-290 


6 


8 


15 


29 


15 


14 


290-300 


6 


9 


14 


29 


13 


15 


0-300 


137 


203 


341 


681 


341 


343 


(lO)). ((Oil. ((671) 


137 


203 


340 


680 







with separation on the left as on the right, where the distribution clearly conies 
from a single ensemble. Recall that Table 1 shows that the frequency distribution 
for zeros is the same for C(0, 1; s) and C(l,4; s). This makes the strong difference 
in the distributions of the gaps all the more interesting. 
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7. Relation between Distributions of Zeros 

We now investigate the relationship between the distributions of zeros of C(l, Am; s) 
and C(0, 1; s). To do this, we consider, rather than the product of these functions, 
as in A3(2, 2m; s), their quotient 

A,i2,2m;s)^'-^^. (7.1) 

C(0, l;s) 

We will now prove some of the important properties of this function. 

Theorem 7.1. The analytic function A4(2, 2to; s) obeys the functional equation 

A4(2, 2m; s) = A4(2, 2m; 1 - s)J-2m{s). (7.2) 

It has first order poles at s — —(2m — 1), —(2m, — 2), ... , —1 and a first order zero 
at s = 0. Its only possible essential singularity is at infinity. On the critical line, it 
lies in either the first or third quadrants, with its argument being given by 



arg[A4(2, 2m; 1/2 + it)] = 02m,c(t) 







(7.3) 



As (T — )• oo for any t, the argument of A4 tends to zero exponentially, while as 
a — ?• —00 for any t,the argument of A^ tends to zero algebraically. 

Proof. The functional equation (|7.2p follows readily from the functional equations 
p.ip for C(l,4m; s) and C(0, 1; s). The poles on the negative real axis arise from 
the poles of J^2m{s), while the zero at s = 1 arises from the first order pole there of 
C(0, 1; s). The function C(0, 1; s) has an infinite set of zeros on the critical line which 
it inherits from ((s) (Titchmarsh & Heath-Brown, 1987). This will make the point 
at infinity an essential singularity unless all but a finite number of these coincide 
with zeros of C(l, 4m; s). The phase of A4(2, 2m; s) on the critical line follows from 
(17^ when 1 - s = s. 

For a positive and not small, we may expand C(l, Am; s) and C(0, 1; s) by direct 
summation: 

C(0,l;.) = 4(^l + l + l+2l + ...), (7.4) 

and 

C(l Am- s) = a(i + f-lV"! + - + g [cos(4mcos-Hl/V5))+cos(4mcos-i(2/V5))] 
V ' ' '' I ^ V '' 2^* 4* S'' 



Combining (17.41) and (I7.5p . we find for m = 1, 2, . . 



(7.5) 



A4(2,4m-2;s)==l--^ + 0(^), (7.6) 



and 



A4(2 4m- s) = l+ ['^os(8"^cos-i(l/V5)) +cos(8mcos-i(2/V5)) - 2] ^ ^^^^ 

T7.7) 
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Thus, in either case, A4 tends to unity in exponential fashion as a increases pos- 
itively. The difference from unity in the case of order 4m — 2 tends to zero as 
exp(— crlog(2)), and in the case of order 4m as exp(— crlog(5)), so that the phase of 
A4 tends to zero exponentially. 

For a negative and not small, we may take the phase of A4(2, 4m; 1 — s) to be 
zero, and write 

Arri^t 
arg[A4(2,4m;s)] ~ arg[J-2„(s)] = 25R[02™(s)] = J^^YJ^jT^^ + Oi^t')- (7-8) 

This quantity tends to zero algebraically as \s\ increases. Its maximum value is 
4?TT,^/(l/2 — (t), which occurs when t ~ 1/2 — a. D 

If necessary, it is not hard to construct uniform bounds on the remainder terms 
in series like that of (|7.5I) for C(l,4m,;s) in the region of absolute convergence 
(T > 1. Indeed, suppose we break the series there into a part J^Ar(l,4m;s) with 
VH+pI — ^i ^H-d a remainder term TIn{^, 4m; s) from ^/pi+pf > ^■ 

|C(l,4m;s)|< | J-^(l, 4m; s)| + |7^^^(1, 4m; s)|, (7.9) 

and I C08{4:m6p-^^p^)\ < 1 for every m, pi, p2, so that 

|-Fjv(l,4m;s)| < |^jv(l, 0; s)|, |7ejv(l,4m; s)| < |7^w(l, 0; s)|. (7.10) 

In turn, 

|J-Ar(l,0;s)| < |J-Ar(l,0;a)|, |7ejv(l, 0; s)| < |7^w(l, 0; a)|, (7.11) 

so bounding uniformly the truncated series, the remainders and the C(l,4m;s). 

Theorem 7.2. The only lines of constant phase of the A4(2,2m;s) which can 
attain a = 00 are equally spaced, and have interspersed lines of constant modulus. 
All such lines of constant phase either reach the critical line in the asymptotic region 
oft at a pole or a zero of A4(2, 2m; s), or curve up and asymptote to each other at 
infinity. 



Proof. From (|7.6p . we have 



-it log 2 



A4(2,4m-2;a + zt) = l-^^^+0(^^), (7.12) 

and so the leading order estimate gives the lines of phase zero for A4(2, 4m— 2; a+it) 
as occurring at i = n7r/log2, for n — 0,1,2,.... Halfway between these lines of 
phase zero are the lines on which the leading order estimate gives 9A4(2,4m — 
2;a + it)/dt — 0: these are lines of constant modulus, and in fact correspond to 
|A4(2, 4m,— 2; a + it)\ = 1, for a — > 00. The same argument applies to A4(2, 4m; a + 
it) , with log 5 replacing log 2 in the estimate for asymptotic placement of lines of 
phase zero and amplitude unity. 

Next, consider where these lines of constant phase and amplitude can go as a 
decreases towards 1/2. Along the line of unit amplitude, as a decreases curves of 
constant phase emanate from the line at right angles and head above and below 
it towards smaller values of a. We see that the phase of A4(2,2m;s) must vary 
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monotonically along the line of unit amplitude; were it to be otherwise, lines of 
constant phase would form closed loops about the line, constraining the phase to 
be everywhere constant within the closed loop- a contradiction. On either side of 
each line of phase zero, we have one line of unit amplitude whose phase increases 
monotonically as a decreases, and another line of unit amplitude along which the 
phase decreases monotonically. A similar argument applies to the lines of phase 
zero, along which the amplitude of A4(2, 2m;s) either increases monotonically or 
decreases monotonically as a decreases. The phase of A4(2,2to;s) on the lines of 
unit amplitude moves away monotonically from zero, so that these lines can never 
intersect lines of phase zero in the finite part of the plane (although they could 
asymptote towards such lines at s — oo). They cannot return to cr = oo, since their 
phase cannot tend back towards zero. Thus, the lines of constant amplitude and 
phase coming from cr = oo must all cut the critical line, or they could from a certain 
point on all commence to curve up towards i = oo. 

Those lines of constant phase zero which reach the critical line in the region 
where the phase of J^2m(s) is accurately constrained by the asymptotic estimate 
p.l4p can only cut the line at a pole or zero of A4(2, 2m; s). D 

We now give one of the two principal results of this section. 

Theorem 7.3. Suppose C{0,1] s) obeys the Riemann hypothesis. Then C{l,^m;s) 
obeys the Riemann hypothesis for any positive integer m. Conversely, if C(l,'im\s) 
obeys the Riemann hypothesis for some m, then C(0, 1; ) obeys the Riemann hypoth- 
esis. 

Proof. Consider lines Li and L2 vtv t > Q along which the phase of A4(2, 2m; s) for 
a given m is zero. Join these lines with two lines to the right of the critical line 
along which a is constant. Then the change of argument of A4(2,2?7i;s) around 
the closed contour C so formed is zero, so by the Argument Principle the number 
of poles inside the contour equals the number of zeros. Each pole is formed by a 
zero of C(0, 1; s), so if there are no such zeros within C there can be no zeros of 
caZC(l, 4r7T,; s) within C. We can repeat this procedure for all m. 

Conversely, each zero of A4(2, 2to; s) is formed by a zero of caZC(l, 4m; s), so if 
there are no such zeros for some m, then there can be no poles, and hence no zeros 
of C(0, 1; s). These arguments prove the theorem in the region to the right of the 
critical line lying between lines of zero phase of A4(2, 2m\ s), with the result to the 
left of the critical line then guaranteed by the functional equation (|7.2p . 

To complete the proof we need to show that any point in the region a > 1/2, 
i > is enclosed between lines of phase zero of A4(2, 2to; s) coming from a — 00. 
We note that i = is one such line, and that for any a > the infinite number of 
such constant phase lines cannot cluster into a finite interval of t, since that would 
indicate an essential singularity of A4(2, 2m; s) for that a. D 

Corollary 7.4. Assuming the Riemann hypothesis applies toC{\^Am] s) orC(0, 1; s), 
apart from a small number of exceptions, lines of phase equal to or-n o/A4(2, 2to; s) 
which leave the critical line going into a < 1/2 cannot reach the asymptotic region 
in which the phase estimate ^7.8}^ is accurate. Instead, they must loop back to cut 
the critical line. 
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Proof. If t is sufBciently large to ensure , _^y^.,* , ^^ < tt, then no line with phase 
TT can enter the asymptotic region. Lines with phase zero are also excluded. Hence, 
lines leaving a — 1/2 with t large and going left which do not return to the critical 
line must either curve upwards and go to t = oo or curve downwards and cut the 
line t = 0. If the first alternative applies, then the image of this line under s — > 1 — s 
in the functional equation is a line of variable phase running upwards in a > 1/2, 
which must cut lines of phase zero running left to intersect the critical line. Such 
intersection points would have to be zeros or poles of A4(2, 2to; s), contradicting our 
assumption. The second alternative means the line of constant phase must intersect 
the axis i = at one of the 2m — 1 poles on the axis. Thus, the exceptional cases 
are limited to the region t < Am? and to lines passing through the 2m — 1 poles. 

In fact, it may be numerically verified in particular cases whether such ex- 
ceptions do in fact occur. They do not for the three cases we have investigated 
(m = l,2,3). D 

We are now in a position to prove the second of the main results of this section. 

Theorem 7.5. Assuming the Riemann hypothesis applies toC{l, 4m; s) orC(0, 1; s), 
then given any two lines of phase zero o/ A4(2, 2to; s) running from a — oo and 
intersecting the critical line, the number of zeros and poles of A4(2, 2m; s) counted 
according to multiplicity and lying properly between the lines must be the same. 

Proof. We consider a contour composed of the two lines of phase zero, the segment 
between them on the critical line and a segment between them in the region cr >> 1. 
The total phase change around this contour is strictly zero, since the region a >> 1 
has the phase of A4(2,2m;s) constrained: — tt << arg[A4{2,2m; s)] << tt. More 
particularly, if P^ = (1/2, i„) lies at the upper end of the segment on the critical 
line, and Pi — {1/2, ti) at the lower end, the total phase change between a point 
approaching P„ on the contour from the right and a point leaving Pi going right is 
zero. This phase change is made up of contributions from the changes of phase at 
the zero or pole Pu, from the zero or pole Pi, from the N^ zeros and Np poles on 
the critical line between P„ and Pi, and from the phase change between the zeros 
and poles. In this list, the first change is 4i2m,c{tu), the phase on the critical line 
just below P„. (We could also have a phase (j)2m,c{tu) — ti", but it will be easily seen 
that in this alternative case the argument which follows will arrive at exactly the 
same conclusion.) The second change is —4>2m.c{ti), where (f'2m,c{ti) is the phase 
just above t;. Giving zero n a multiplicity z„, and pole n a multiplicity p„, the 
phase change at the former is —nzn and the latter 7rp„. The phase change between 
zeros and poles is 4'2m,c{ti) — 4>2m.c{tu)- Hence, the phase constraint is 

(f>2m,c{tu) - <t>2m,c{tl) - T^\^ ^" ~ X! ^"1 "*" '^2m,c(tO " (f>2m,c{tu) = 0, (7.13) 

n— 1 n— 1 

leading to 

N^ Np 

2z„ = ^p„, (7.14) 

n— 1 n— 1 

as asserted. D 
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Corollary 7.6. // all zeros and poles on the critical line have multiplicity one, the 
numbers of zeros and poles on the critical line between any pair of lines of phase 
zero of A4(2, 2to; s) coming from a = oo are the same. The distribution functions 
for zeros -/Vci,4(2,i) o,rid -/Vco,i(2'^) '^f (^3' ^^S'^ must agree in all terms which go 
to infinity with t. 

Proof. The first assertion is a simple consequence of Theorem 7.5. The second 
follows from Theorem 7.2 and Theorem 7.5: the number of zeros and poles between 
successive zero lines coming from a = oo match for all such pairs of lines, and 
there are only a finite number of exceptional poles and zeros which may disturb the 
equality between numbers of zeros and poles. If lines of phase zero curve up and 
asymptote to the point at infinity, we apply Theorem 7.5 to the leftmost of these 
and the last line of phase zero cutting the critical line. D 

We now show in Fig. [T3] and [TJ] some examples of morphologies of lines of 
constant phase of A4(2, 2; s) near the critical line. Fig. [T51 gives at top left a diagram 
in which there are three lines of phase zero coming from a = oo. The first two give a 
loop going right which intersects the critical line at a pole, below, and a zero above. 
The third again intersects the critical line at a pole, and ends what we will describe 
as a "cell" . By this we mean what is a repeat unit in the loose sense; above the third 
line we start another loop going right. Fig.[l3]shows a view at a structure for higher 
t in the upper right graph, in which the resolution of contour lines is insufhcient to 
give a correct impression of the relationship between the lines of constant phase. 
The magnified views below show at left that the apparent intersection of lines of 
phase and tt is in fact an avoided crossing, where the lines behave in hyperbolic 
fashion. In the graph below right, we see a small phase loop going left away from 
a = 1/2, with an even smaller line element of phase occurring in a > 1/2. Note 
that, denoting poles and zeros by prefixes P and Z, the values of t at which zeros 
and poles occur in the upper graph are: 

P45.6, Z45.9, Z46.9, Z47.71, P47.74, P48.0, Z49.2, P49.72, P49.77. 
This data then shows a succession of three zeros of C(0, 1; s) uninterrupted by zeros 
of C(l,4;s). We have encountered cases where the opposite situation occurs, but 
no sequences of four successive zeros or poles. 

The leftmost graph in Fig. [14] shows the hyperbolic distribution of lines of con- 
stant phase occurring around a zero of A'(2, 2; s). The discussion following Theorem 
5.6 may be extended from A3(2, 2m; s) to A4(2, 2m; s), with the result that between 
successive zeros of A4(2, 2m; s) uninterrupted by a pole there will be a hyperbolic 
point to the right of the critical line. Between successive poles of A4(2, 2m; s) un- 
interrupted by a zero there will be a hyperbolic point to the left of the critical line. 
The leftmost graph shows the case of two poles, with a hyperbolic point in a < 0.5. 
The rightmost graph in Fig. [HI shows the morphology of lines of constant phase for 
A4(2,4;s), for which the density of lines of phase zero coming from a = oo is, as 
expected, greater than for A4(2, 2; s). 

Table 2 gives the numbers of zeros and poles in successive cells, for the first 
twenty- two cells of A4(2, 2; s). These evidently satisfy the requirements of Theorem 
7.4. Note the tendency of the numbers of zeros and poles to increase with cell 
number; this is a consequence of the cell length tending to be roughly constant, while 
the density of poles and zeros increases logarithmically, in keeping with equation 
(16.71) . For the fourth column, the mean length is 9.01 and the standard deviation 
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Figure 13. Contours of constant phase for A4(2, 2; 

±7r in green. 



with phase zero in blue, and phase 



is 0.90. Note that the corresponding separation at cr = oo is 27r/log(2), or around 
9.06, so that there is no significant change in average spacing with variation of a. 
For A4(2, 4; s), Table 3 gives data for the first twenty-four cells, with the mean cell 
length being 3.92 and its standard deviation 0.53. The separation at infinity here 
is 27r/log(5), or around 3.90. 

Corollary 7.7. // all zeros and poles on the critical line have multiplicity one, 
on any closed contour (including those closed at a ^- oo) there must be an equal 
number of zeros and poles and they must alternate, i.e. following the contour in a 
clockwise direction, any pole will be followed by a zero and zero will be followed by 
a pole. Every cell begins with a pole on the critical line. 

Proof. We follow the line of zero phase coming in from a = oo which forms the 
start of the cell. It must intersect the critical line, in accordance with Corollary 7.2 
and meet there a line of constant n phase going towards increasing t, which is also 
a branch cut discontinuity. The phase 4>2m,c{t) is decreasing as t increases along 
the critical line and above the intersection point (call it ii ) there is a change of tt 



Article submitted to Royal Society 



30 



R.C. McPhedran and others 




Figure 14. Contours of constant phase for A4(2, 2; s),(left) near a hyperbolic point, 
(right) with phase zero in blue, and phase ±7r in green. 

to the phase on the critical hne, forcing the phase to drop to 4'2m,c — t^, which lies 
between (— tt, 0). As phase is increasing in a clockwise direction around the point 
of intersection, that point must be a pole (i.e., every cell begins with a pole). 

We continue to follow the branch cut as it loops back to the critical line in 
accordance with Corollary 7.4. The phase below the next intersection point t2 is 
still negative as <t>2m,c{t) decreases along the critical line. The change of tt to the 
phase at ^2 niust be a positive change, such that the phase above the point of 
intersection must be 4'2m,c{i2) > . The phase is increasing in an anti clockwise 
direction, which implies that the intersection point ^2 is a zero. This argument 
is continued for successive poles and zeros until we reach the end of the cell in 
question. D 

Remarks: 

The following consequences are clear from the arguments in Corollary 7.7: 

• Any closed contour intersecting the critical line at a finite set of points will 
have maximum and minimum intersections of opposite type i.e a zero and a 
pole. 

• If there are two consecutive zeros or poles on the critical line they cannot lie 
on the same closed contour of zero phase. 

• Any closed loop of zero phase going from cr = 1/2 into a < ^ must have a 
zero at its lowest intersection with a = 1/2. 

The research of R.McP. on this project was supported by the Australian Research Council 
Discovery Grants Scheme. 
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Table 2. Numbers of zeros and poles, interval m t, and starting point on the critical line, 
for the first twenty-two cells lying between lines of phase zero coming from a — oo, for 

A4(2,2;s). 
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# Zeros 


#Poles 
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Lowest pole 
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3 


3 


8.11 
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2 


5 
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14.13 


3 


7 
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9.66 


23.28 


4 


6 


6 


7.98 


32.94 


5 


8 


8 


8.85 


40.92 


6 


9 


9 


10.65 


49.77 


7 


8 


8 


7.95 


60.42 


8 


9 


9 


8.77 


68.37 


9 


10 


10 


10.29 


77.14 


10 


10 


10 


8.71 


87.43 


11 


9 


9 


8.19 


96.14 


12 


11 


11 


9.99 


104.33 


13 


10 


10 


8.63 


114.32 


14 


10 


10 


8.14 


122.95 


15 


12 


12 


10.16 


131.09 


16 


11 


11 


9.05 


141.25 


17 


10 


10 


8.41 


150.30 


18 


11 


11 


8.47 


158.71 


19 


13 


13 


10.52 


167.18 


20 


11 


11 


7.90 


177.70 


21 


11 


11 


8.63 


185.60 


22 


13 


13 


9.99 


194.23 
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Table 3. Numbers of zeros and poles, interval in t, and starting point on the critical line, 
for the first twenty-four cells lying between lines of phase zero coming from a = co, for 

A4(2,4;s). 



Cell# 


# Zeros 


#Poles 


Length 


Lowest pole 


1 


1 


1 


4.22 


6.02 


2 


2 


2 


3.89 


10.24 


3 


2 


2 


4.16 


14.13 


4 


2 


2 


3.16 


18.29 


5 


3 


3 


4.28 


21.45 


6 


2 


2 


3.93 


25.73 


7 


3 


3 


3.28 


29.66 


8 


3 


3 


4.65 


32.94 


9 


3 


3 


3.33 


37.59 


10 


4 


4 


4.68 


40.92 


11 


4 


4 


4.17 


45.60 


12 


3 


3 


3.20 


49.77 


13 


3 


3 


3.96 


52.97 


14 


4 


4 


3.90 


56.93 


15 


4 


4 


4.28 


60.83 


16 


3 


3 


3.26 


65.11 


17 


4 


4 


3.79 


68.37 


18 


4 


4 


4.54 


72.16 


19 


4 


4 


3.51 


76.70 


20 


5 


5 


4.53 


80.21 


21 


3 


3 


2.89 


84.74 


22 


4 


4 


4.61 


87.63 


23 


5 


5 


3.90 


92.24 


24 


4 


4 


4.00 


96.14 



Appendix A. Supplementary Notes on Systems of Angular 

Sums 

We give here results linking trigonometric lattice sums of order up to ten, which 
show that they may be generated from three independent sums, C(0, 1; s), C(l, 4; s), 
and C(l, 8; s). We also give expressions in terms of sums of Macdonald functions of 
the Kober-type from which these independent sums may be calculated. 
A basic result we use relies on the symmetry of the square lattice: 

C(l, 4m - 2; s) = ^ ' i!!^ 7 J\^'^^^ = ^ = '^(^^ 2™ - 1; s) - 5(2, 2m - 1; s), 



{pI+pIY 



(Al) 



from which we obtain 



C(2, 2to - 1; s) = 5(2, 2m - 1; s) = -C(0, 1; s). 



(A 2) 



for m any positive integer. We can also expand cos(4?ti— 2)6'pj^^p2 — T4m-2{cos9p-^^p2) 
in terms of powers of cos 9pi,p2 , using the expressions for the Chebyshev polynomials 
in Table 22.3 of Abramowitz & Stegun (1972). This enables us to express C(4m — 
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2, 1; s) in terms of C(4n, 1; s) with An < Am — 2, and n being a positive integer. In 
this way, we can inductively arrive at the results below. 

(a) Order 

The only sum of this type is C(0, 1; s), given by equation (J2.2I) . It is given by a 
modified form of (j2.3p . since there are contributions from both axes pi — 0, p2 — 
rather than just p2 — 0: 

C(0,l;.) = 2C(2.)+ '^^^;~^/') c(2.-l) 

T{s) 

_ ^ oo °o / \ s-1/2 

■ E E (?) K,-iM2np,P2). (A3) 



r(s) ^^ ^^ \pi 

(6) Order 2 



The single sum of order 2 is 



2V^r(. + i/2)c(2,s-i) , 87r- ^^^p,y-v^ 
= ri^TT) + rxTTT) -^ ^[vj ^i^2'^^^+V2(2-piP2). 

(A 4) 

(c) Order 4 
We generate this system from C(4, 1; s), obtained from (j2.3p . In terms of this, 

C(l, 4; s) = 8C(4, 1; s) - 3C(0, 1; s), (A 5) 

and 

C(2, 2; s) = 4C(4, 1; s) - C(0, 1; s), 5(2, 2; s) = -4C(4, 1; s) + 2C(0, 1; s). (A 6) 

(rf) Order 6 
The result of expanding (IA1|) for rri = 2 is 



C(6,l;s) = ^C(4,l;s)-ic(0,l;s), (A 7) 

while from (|A2I) 

C(2,3;s)=5(2,3;.s) = ic(0,l;s). (A 8) 

Also, 

v^ , p\pI \p / pIp\ ^ 1 v^ ^ PiPJ ..q^ 

/- (p2+p2)s+3 Z. (p?+Pi)^+' 2/- (P?+Pi)^+'' ^ ^ 

(Pl:P2) (Pl,P2) (Pl,P2) 
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Figure 15. The modulus of C(l, 4; s) (red) and C(l, 8; s) (blue) as a function of 
s = l/2 + it, for t e [0, 20] (left) and t G [20, 40] (right). 



SO 



E 

(Pl,P2) 



, cos^(6>p,,pJsin^(6lp,,pJ ^ 1^ ^^ ^ 

(p2+p2)s 8 ^ ' ' ^" 



(A 10) 



Note that all sums mentioned so far can be generated from just two, say C(0, 1; s) 
and C(l, 4; s). C(6, 1; s) is given by (^3]) . 



(e) Orrfer ,S 
The new sum we use here is C(8, 1; s), obtained from (12.31) . In terms of this, 

C(l, 8; s) = 128C(8, 1; s) - 224C(4, 1; s) + 49C(0, 1; s), (A 11) 

C(2, 4; s) = 64C(8, 1; s) - 112C(4, 1; s) + 25C(0, 1; s). 



and 



5(2, 4; s) = -64C(8, 1; s) + 112C(4, 1; s) - 24C(0, 1; s) 
Two other sums in this system are 



E' 



pfpi 



(P1.P2) 



(PI+P2 



,2\s+4 



-C(8,l;s) + ^C(4,l;s)-ic(0,l;s), 



and 



E' 

(Pl:P2) 



4 4 
P1P2 



(p2^p2).+4 = C(8, 1; s) - 2C(4, 1; s) + -C(0, 1; s). 



(A 12) 
(A 13) 

(A 14) 
(A 15) 



Using the result (|A 11[) . wc have calculated curves showing the modulus of 
C(l,8;s) as a function of s = 1/2 + it on the critical line. These are given in 
Figs. [TMTBl and the distributions of zeros they show are given in Table [TJ 



(/) Recurrence Relations 

We can generalize the above procedure by establishing recurrence relations for 
the trigonometric sums. We consider 



E' 

(Pl,P2) 



Pf" 



(pI+pI) 



,2\s+n 



C(2n,l;s) 



E 



1 



(Pl,P2)#(0,0) 



(pI+pD 



1 



P2 



pI+pI 



(A 16) 
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Figure 16. The modulus of C(l,4;s) (red) and C(l,8;s) (blue) as a function of 
s = 1/2 + it, for t € [40, 60] (left) and t £ [60, 80] (right). 



Expanding using the Binomial Theorem, we obtain 



C{2n,l;s) ^J2"Q{-iy 0(21,1; s). 



(A 17) 



1=0 



If n is even, (IA17P gives an identity: 



2n-l 



Cii-iy 0(21,1; s)^0, 



(A 18) 



while for n odd we obtain an expression for C(4n — 2, 1; s) in terms of lower order 
sums: 



_. An — A 

0{An - 2, 1; s) = - 5] ^n-i^^ (_^y ^(3/, 1; s). 



(A 19) 



1=0 



It may be checked that the two relations (lAlSp and (|A 191) give equivalent results 
in the cases given above. 

In Fig.[T7]we give null contours of the real part and imaginary part of A3(2, 4; s), 
for comparison with those of A3 (2, 2; s) in Fig. [3] It will be noted that there are 
now three lines giving the contours on which ImA3(2,4;s) = 0, and these are not 
circles (as was the single exemplar in Fig. [3]). There is a corresponding increase 
in the number of null contours starting and finishing on the real axis, and in the 
value of (T at which the last one reaches the real axis. There are two examples 
in Fig. [T7] of zeros on the critical line of the real part, but not the imaginary part 
(tan ('/>2m,c{l/2 + it) = cxd), and one of a zero of the imaginary part but not the real 
part (cot (02771, c(l/2 + it) ~ 00). The null contours starting at cr = cxd and ending 
at cr = —00 settle down to an asymptotic behaviour for somewhat larger values of 
t than in Fig.[3J but are otherwise similar to the lower order null contours. 



(5) Order 10 
The new sum we use from (|2.3p is C(10, 1; s), which, from (jA 19[) . is 

C(10, 1; s) = ic(0, 1; ,s) - ^C(4, 1; ,s) + ^C(8, 1; s). 



(A 20) 
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Figure 17. Null contours of the real part (red) and imaginary part (blue) of 
A3(2, 4; cr + if), with a G [-4.5, 11.5], and t £ [0.1, 20]. 



Other 


sums 


in 


the 


system are 














(Pl>P2) 


{vl 


pIpI 

+Piy 


+5 


and 








\ ^ ' 




pIp\ 





C(8,l;s)-C(10,l;s), 



- E '- 



4 4 
PiP2 



(Pl,P2) (Pl,P2) 



(A 21) 



(A 22) 



Appendix B. Derivation of the Functional Equation for 

C(l,4m;s) 

We start with the equation (30) from McPhedran ei al (2004). This equation uses 
the Poisson summation formula to derive a connection between a sum over a direct 
lattice of points Rj, = Rpi.pa = d{pi,p2) with polar coordinates {Rp^4>p), and a 
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corresponding sum over the reciprocal lattice, where the lattice points are labelled 
K?j = 'Khi,h2- Considering the case of a square lattice with period d, the reciprocal 
lattice points are Kh = (27r/d)(/ii, /i2)- The sum in the direct lattice incorporates 
a phase term of the Bloch type, with wave vector ko, and the result of the Poisson 
formula is 



2^^-^r(i+.) y: ' 



2 ' ' ^^—' i?2s 

{Pl;P2) P 

= ^r(- + l-.) ^' |. (Bl) 

Here the sum over the reciprocal lattice in fact runs over a set of displaced vectors: 
Qh = ko + K,, = (^ + fcor., ^ + k^y) = [Qh,eh). (B 2) 



where in (IB 2\ the second and third expressions are in rectangular and polar coor- 
dinates. Note that in (jB II) and (IB 21) . pi,p2 and hi, ft,2 run over all integer values, 
and 9o gives the direction of kp. 

We express the relation (IB II) in non-dimensionalized form, taking out a factor 
d^* on the left-hand side, and a factor [2ti / d)'^^'^^ on the right-hand side. We put 
ko = (27r/d)/t0 5 a-nd obtain 

r(, + »)E " (,;,,.)/ -iV-r(, + i-.) 

(Pl,P2) 

If we now let ko — ^ 0, set / — Am and take Re(s) > 1, we obtain the desired result: 
r(2m + s)^, , T(2m + l- s)^, , , , 
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